Electrons and holes in graphene under electromagnetic waves: gap appearance 

and non-linear effects 
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We find the exact solution of graphene's carriers under electromagnetic radiation. To obtain such 
solution, we combine Floquet theory with a trial solution. Then the energy spectrum is obtained 
without using any approximation. Using such results, we prove that the energy spectrum presents 
a gap opening which depends on the radiation frequency and electric wave intensity, whereas the 
current shows a strongly non-linear behaviour. 
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Graphene, a two dimensional allotrope of carbon, 
has been the centre of a lot of research since its 
experimental discovery four years ago It has 
amazing properties (HQ] HI- For instance, electrons 
in graphene behave as massless relativistic fermions 
. Such property is a consequence of its bi- 
partite crystal structure 0, in which a conical dis- 
persion relation appears near the K, K' points of 
the first Brioullin zoneQ. Among other properties 
one can cite the high mobility that remains higher 
even at high electric-field induced concentration and 
translates into ballistic transport on a submicron 
scale at 300°K. Graphene is therefore a promis- 
ing material for building electronic devices, but there 
are some obstacles to overcome. One is the trans- 
mission probability of electrons in graphene, which 
can be unity irrespective of the height and width of 
a given potential barrier. As a result, conductivity 
cannot be changed by an external gate voltage, a fea- 
ture required to build a field-effect transistor (FET) , 
although a quantum dot can be used 10]. In a pre- 
vious paper, we have shown that a possible way to 
induce a pseudogap around the Fermi energy con- 
sists in doping graphene (Tlj . On the other hand, 
much efforts have been devoted to understand the 
electrodynamic properties of graphene as well as its 
frequency dependant conductivity , [l3[ , [l4[ , [l5[ . 
In this Rapid Communication we solve, the prob- 
lem of graphene's electron behavior in the presence 
of an electromagnetic plane wave. As a result, in 
the limit of long wave length of the electromagnetic 
field, we are able to find a gap opening. This is like 
if electrons in graphene acquire an effective mass un- 
der electromagnetic radiation. We also calculate the 
current and show that there is in fact a strongly 
non-linear electromagnetic response, as was claimed 



FIG. 1: (Color online). Graphene lattice. Atoms in the 
A sublattice are shown with different color than those in 
the B sublattice. The vector k in which the electromag- 
netic field propagates and the electric field vector E are 
shown. 



before by S. A. Mikhailov[l6j using a semi-classical 
approximation. 

Considering an electron in a graphene lattice sub- 
ject to an electromagnetic plane wave as shown in 
Fig. 1. The plane wave propagates along the two- 
dimensional space where electrons move. In this 
Rapid Communication we take k = (0,fc). The 
generalization to any direction of k is straightfor- 
ward. It has been proved that a single-particle 
Dirac Hamiltonian can be used as very good ap- 
proximation to describe charge carriers dynamics 
in graphene fl7l|. For wave vectors close to the K 
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point of the first Brioullin zone, the Hamiltonian is 

ta, 

where vp is the Fermi velocity vp « c/300, 7r = p — 
eA/c with p being the electron momentum operator 
and A is the vector potential of the applied elec- 
tromagnetic field, given by A = (— a cos(ky — tut), 0) 
, where Eq is the amplitude of the electric field and 
to is the frequency of the wave. E is taken as a con- 
stant since screening effects are weak in graphene 
[l~9l |. For the valley K' the signs before ir y are the 
opposite [l8[ . The dynamics is governed by 
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+ [£ 2 cos 2 4>—£,vpha z k sin (f>—ifuv^a x sin 0] = 0. 
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where we have defined the phase <f> of the electro- 
magnetic wave as <fr = ky — Lot. The parameter £ 
is defined £ = e ^J F , and cr^ is the set of Pauli 
matrices. To solve this equation, we follow Volkov 's 
suggestion to use a trial function that depends upon 
the phase 4> of the wave with the following form, 



H(x,y,t)^(x,y,t) =ih 



dt 



(2) 



where 



is a two component spinor. Here A and B stand for 
each sublattice index of the bipartite graphene lat- 
tice Q-To find the eigenstates and eigenenergies we 
adapt a method developed by D. M. Volkov in 1935 
[20| to study the movement of relativistic particles 
under an electromagnetic field. Volkov considered 
the complete four component bispinor and Dirac ma- 
trices. We use instead only Pauli matrices, and two 
component spinors. First we write the equations of 
motion for each component of the spinor 



y, t) = e lp * x/h+lp y y/h - iet/h ¥{(j>), (8) 



where e = vp^jp x + P 2 and F(<p) is a spinor. In- 
serting Eq. ijHJ) into Eq. ( [TJ) , it yields a equation for 
F(</>). The resulting differential equation is, 



-h 2 (v F k 2 -u 2 ) 



-2iherj — jj^-+[—2£vpp x cos (f>—^vpha z ki 
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+ £ 2 cos 2 <j>- ihw^a x sin <f>]F(<f>) = 0, (9) 

and r\ — clo — vpkp y . The solution to previous equa- 
tion has the form F(0) = exp(i<j)B /2A)z(<p), where 
z((j)) is the solution of 



v F (it x - nTy)W B {x,y,t) = in — , (3) 



0t 



vf{tt x + iTT y )^ A (x, y, t) = iK 



&SF B (x,y,t) 
dt 



(4) 



Considering the magnetic field as B, the commuta- 
tion rules for tt x and tt v are in this case, 



[7ti,Ttj] = — EijfcBfc, i,j = x,y (5) 



z'\4>) + [{C{4>)/A) - B 2 /(2A 2 )]z(0) = 0, 

A = -h 2 (v 2 F k 2 - lo 2 ), B = 2hi] and 
C(<p) — —2£vfp x cos <f> — ^vpha 2 fcsin <f> + £ 2 cos 2 <\> — 
ihw^a x sin </>. On the other hand this equation can 
also be solved if we consider the long wave limit 
(A — * oo) and if it is neglected the term with sec- 
ond derivative upon <j> to give, 



P(t)=exp[G(t)]u, (10) 
where u is a two component spinor and 



— , it x db ifc v 
dt 1 v 



eE a 

= sm(ky — Lot), (6) 

c 



and using that k^A^ = 0, we find the following equa- 
tion of motion for the spinor: 



G(t) = -—7-t - — — sinu)tp x 
4ne neuj 

+ —$ — sm2tot — a x cos tot, (11) 
Sneoj 2e 

The complete solution to our problem is: 
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and $>(x,y,t) is the Floquet mode Eq.f [T5|) .I f the 



^(x, y, t) = exp [ip x x/h + ip y y/h — iet/H] exp [G(t)] u, 

(12) ' 

where u must be taken as a two component spinor 
which satisfies the requirement that when A — > 0, 
^(x,y,t) should be the solution of the free Dirac 
equation to avoid strange solutions. Thus u is given 
by: 



U = 7l( ± 



p = tan(^) (13) 

Px 



where the sign plus stands for electrons and minus 
sign stands for holes. The eigenvalues are e + £ 2 /4e. 
But now we need to take into account the time pe- 
riodicity of the electromagnetic field, which means 
that the solution must be written in the following 
way P, 



* = exp(-ixt/h)${x,y,t), 



(14) 



where y, t) is periodic in time, i .e. *&(x, y, t) = 
<&(x,y,t + T), and \ is a re& l parameter, being 
unique up to multiples of Huj, u) = 2tt/t. In fact, 
such property is very well known within the so-called 
Floquet theory [2l| . developed for time dependent 
fields. This time periodicity of the field leads to the 
formation of bands. Our solution Eq. (fT2"|) satisfies 
the requisite form Eq. (fT4")l . However, to obtain the 
Floquet states plj for this problem, we need to con- 
sider that e + £ /4e is unique up to multiples of hw. 
Thus E — e + £ 2 /4e+ nhuj : where n is an integer 
such that n = 0, ±1,±2, while <fr(x,y,t) is the 
Floquet mode given by, 

&{x, y, i) = exp [intot + ip x x/h + ip y y/h] exp [-F(i)] u. 

(15) 

where F(t) = G(t) + if t/AHe. The Floquet modes 
of the problem must satisfy the equation [2l| . 



H (x, y, t) — ih 



*(x,y,t) = E3>(x,y,t) (16) 



The final cigcncncrgics for electron and holes are 
the following, 



E n (p) = nho ± v F \\p\\ ± 



e 2 Elv F 



Ac 2 u> 2 ||p|| ' 

where n — 0,±1, and the wave-function is, 



(17) 



two smallest terms are neglected in Eq. ( ITT|l. our so- 
lution can be reduced to one found in Ref. [22] • The 
spectrum given by Eq. (|17[) is made of bands, where 
the electromagnetic field bends the linear dispersion 
relationship due to the last term. 



A(p) = 



e 2 Elv F 



P 



£ 

4c 



(19) 



Thus, around the Fermi energy, holes and electron 
bands are separated by a gap of size A(p) = 2A(p F ). 
To estimate the magnitude of the gap, we consider 
electrons near the Fermi energy, thus e = v F \\p\\ « 
e F , where e F is approximately [l3j e F = 86meV. 
For a typical microwave frequency lo = 50GHz with 
an intensity Eq = iV/cm of the electric field, the 
gap size is around A « 0.2meV . Due to the gap 
opening, the particles are not longer massless. The 
mass acquired by the carriers due to the electromag- 
netic field is therefore around 10~ 4 m e . Since this 
is a time-dependant problem, one cannot measure 
the gap directly from the density of states. Instead, 
one can look for jumps in the conductance, using, 
for example, the device proposed in Ref .[24]. Fi- 
nally, the solution can be used to evaluate the cur- 
rent using the collisioness Boltzmann equation [l6|, 
neglecting interband transitions. However, this ap- 
proach turns out to give the same results as con- 
sidering the velocity of the particles as time depen- 
dant while the distribution function [f p /t)] remains 
static [2 3]. Such equivalence is the result of the triv- 
ial statement that in the Boltzmann approach, the 
electric field induces a displacement fo the Fermi 
surfacc[23].In fact, the current obtained below is 
equal to the one obtained using a Boltzmann equa- 
tion semiclassical approach for graphene under an 
electric adiabatic field pjj], which also reproduces 
the intraband Drude conductivity [Hj]. The electric 
current is given by j(t) = AS^ 1 E P Jp(*)/p(*)> where 
5" 1 is the sample area and the factor 4 comes for the 
spin and valley degeneracy. j p (i) denotes the con- 
tribution to the current of particles with momentum 
p at time t, while f p {t) is the distribution function 
in phase space at temperature T. Let us first cal- 
culate the \i component of the current vector, given 
by j^{t) = ev F ^* iP a^ niP . We use Eq. (18]) to 
obtain the components of j in the x and y directions, 



*„ iP (x, y, t) = exp(-iE n (p)t/h)$(x, y, t), (18) 



jx,p{t) — ev F sinh ^— cosaji^ +cosh ^-coswt^ cos <p, 

(20) 

and jyjp — ev F simp . In j x we already observe a 
very important non-linear behavior. In fact, j x ,p(t) 
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can be written as a combination of harmonics by 
using a Fourier series development of the hyperbolic 
functions [2^1. 



jx, P (t) = ev F ^ J 2s+i (7 ) cos((2s + l)uxb) 

s=Q ^ 



a 



where Qo 
obtain, 



(1 + 



£/ef = zEq/ '(cajpp). For s — we 



ev F I Jo ( - 



2^ J 2s ^-^ cos(2sa;i)^ cos Lp, 



. (21) 

where J s (^/e) is a Bessel function. Our result is sim- 
ilar to the first terms of the current obtained using a 
semi-classical approximation [l6| , except for the last 
term which eventually cancels out in the thermody- 
namical limit. Now we include such limit by using 
f P (t). In this case, we are dealing with quasiparti- 
cles with an effective dispersion relation. Using Eq. 
([2"Tj) . and by summing over the phase space we get, 

A 00 



00 00 



■hs+i [ — ) n(p)dp x dp y , 



where n(p) = [1 + exp(E n (p) — ep)/ksT}~ 1 is the 
occupation factor. For ep 3> fc^T, n(p) can be re- 
placed by a step function. Using Eq. (jTTJ) , v F p F — 
e/(l + VI - (£A/) 2 )/2 for n = 0, A(s) is, 

1 

A(s) = 2^a 2 y J 2s+1 .tc?x, (23) 



A(0) w 7rp|aQo 1 



Qo 

a 



1 
576 



(24) 

and for s > 0, the integral can be approximated as 
A(s) w 2ttp| (Q /2) 3 r(s - |)/ar(s + |). The first 
terms of the current are, 



j x (t) w en e v F Q 



(1 - jQg) cos(^) 
-^Qo cos(3wt) 



where n e is the density of electrons n e = pp/irh 2 . 
In conclusion, we solved the Dirac equation for 
graphene charge carriers under an electromagnetic 
field. 
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